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Variations of polarisation in external electrostatic fields 
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Abstract. - It is nowadays a quite diffuse idea that variations of polarisation in condensed matter 
theory are related to a "Berry phase". The derivation of the latter geometric phase is correct 
only if the restrictive periodic gauge [1] is used for its derivation, aside giving rise to a "quantum 
of polarisation" that is of difficult theoretical interpretation. Its derivation has not been demon- 
strated in the general case of an external homogeneous electric field interacting with the electronic 
field. In the present paper we give a brief derivation of the polarisation differences in a general 
manner than in [1] allowing for a general boundary conditions freedom showing that is possible 
to define it even in the presence of an external electric field interacting with the electronic field 
and writing equations in a way directly implementable in modern first principles codes. 
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[ Introduction. — A revolution of old concepts of clas- 
sical electrostatics has been achieved by the work of sev- 
eral authors [1, 2], where we cite in particular that of 
Resta [3,4]. There the concept of variations of polarisation 
AP as a geometrical quantum phase has been introduced, 
and only the definition of the polarisation of the electronic 
field in crystalline materials in absence of external electric 
field has been considered. Following a different approach, 
we shall define the polarisation differences in a quite gen- 
eral manner allowing for its calculation when it is the case 
of an external homogeneous electrostatic field interacting 
with the electronic system. Also, in our derivation we 
still consider changes of AP which occurs upon making a 
change in the Hamiltonian (Kohn-Sham Hamiltonian if we 
consider a density functional theory context) of a physical 
system writing equations in a way directly implementable 
in modern first principles codes. 

Polarisation differences. — As it is already well 
known [2] the electric dipole of a macroscopic but finite 
piece of matter is not a bulk property as the latter de- 
pends upon truncation and shape of the chosen sample. 
We agree that only variations of the latter are experimen- 
tally accessible where it is supposed that the state of the 
physical system is specified by a certain parameter A. We 
consider as a formal definition of the electronic dipole the 
following expression: 
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total electric charge due to the electronic and ionic con- 
tributions. Making use of the Born-Oppenheimer approx- 
imation we can divide the electronic and ionic contribu- 
tions in expression (1) and express the electronic charge 
contribution as follows: 



Pel = e 



dk 
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where |4'n,k > are the eigenstates of the Hamiltonian 
(dependent on A) and /„ are the occupation factors. We 
always assume eigenstates of the Hamiltonian of being in 
the Bloch form: 



|*n,k > = 



k > 
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Within this gauge representation it is naturally appar- 
ent the parametric dependence of electronic eigenstates 
on the wave vector k. Here quasi periodic [5] boundary 
conditions are concerned, where eigenstates |^'„.k > are 
not periodic while |itn,k > are periodic on a lattice and 
orthonormal in the unit cell but not periodic in the recip- 
rocal space. Within this boundary conditions the eigen- 
states arc differentiable and wave vectors k are continuous 
and can be treated at the same foot of the r coordinate. 
On the other hand our derivation of the polarisation dif- 
ferences can equally be brought forward assuming the pe- 
riodicity of |u„_k > in reciprocal space. Let us assume the 
Hamiltonian of the system being of the form: 



where 17 is the volume of the system and p\{r) is the 



F(A) =T + y(Eo,A) 
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where Eq is an external homogeneous electrostatic field. 
We consider by hypothesis that the system is in a 
steady state [6] , and focus our attention only on the elec- 
tronic part of the electronic dipole trying to calculate the 
polarisation derivatives with respect to the parameter A. 
Bearing in mind eq.(2) we can easily determine deriva- 
tives [1,2] of the dipole of the system with respect to A 
in a finite sample and then take the thermodynamic limit. 
The dipole derivative is reported in eq.(5), where e„^k are 
the eigen-energies of the electronic system, c.c. represents 
the complex congiugate, Q the volume of the unit cell, N 
the number of cells of the system, p is the momentum 
operator. 

Up to know our derivation is similar to Resta deriva- 
tion [2] apart from a different choice of the boundary 
conditions for eigenstates \^n,k > said eigcn-amplitudes 
\unM >, and the assumption that the electronic system 
is in a steady state and interacting with an external ho- 
mogeneous electrostatic field. At this point we can recast 
eq.(5) in terms of eigen-amplitudes |un,k > and report its 
expression in eq.(6). There i7k is the effective Hamiltonian 
Hk = e-'^'^He'^". 

Before to proceed further we consider the matrix iden- 
tity already demonstrated in [7] and reported in cq.(7). 
Let us now substitute eq.(7) in cq.(6) and obtain cq.(8). 
At this point we can easily estimate the infinite sum in 
eq.(8) making use of Salem sum rule [8] and recast eq.(8) 
in a more compact and easy calculable form obtaining: 

e d 

^ ''^ ^ 'mi ^ ^ ^ "™.kN^k|un,k > (9) 
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then integrating in A between two different states labelled 
by A = and A = 1 and taking the thermodynamic limit 
we can recast the integral of eq.(9) as follows: 

AP,,=PW-P(°) (10) 

being. 
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that is formally the same result obtained in [1] for the 
case of no electric field interacting with the electronic field. 

Our result it is general and allows for calculations of po- 
larisation differences even in the case of an external elec- 
tric field interacting with the electronic field. Aside, we 
do not need anymore ask for the periodic gauge [1], i.e 
eigen-amplitudes jwn.k > periodic in k. We achieve the re- 
sult that polarisation differences are always calculable by 
eq.(ll) allowing us to make use of a more ficxible bound- 
ary condition. On the other hand our derivation still holds 
even in the case we use a periodic gauge. In turn our re- 
sult allows for the elimination of the so called " quantum of 
polarisation" [1] that is basically due to the very restric- 
tive periodic gauge, attracting attention on the fact that 
the presence of an external electric field breaks the trans- 
lational symmetry in the direction of the field and that 
perhaps the correct boundary condition of the physical 
problem might not be periodic in reciprocal space. 

Summary. — During the analysis of the problem 
started by considering the interaction of the electronic field 
with an external electric field it has been shown how to 
calculate variations of polarisation. We firstly make use of 
the expression derived by Resta for the polarisation deriva- 
tives then we calculate, in a very general manner without 
resorting to the restrictive periodic gauge, the polarisa- 
tion differences. Our result is independent from deriva- 
tion brought forward in ref [1] where it was the case of no 
interaction of the electronic field with an external homo- 
geneous electric field. In that case it was considered the 
polarisation difference as a Berry phase of the system but 
the latter derivation needs of a very restrictive boundary 
condition, the so called periodic gauge. Without the above 
mentioned boundary condition the algorithm proposed in 
ref [1] Jails while our derivation still holds allowing for a 



general gauge freedom and also allowing for the elimina- 
tion of the concept of " quantum of polarisation" , the latter 
being theoretically difhcult to be interpreted. Aside, we 
wrote equations in a way directly implementable in mod- 
ern first principles codes. 
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